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SUMMARY 


This report discusses measurements of soot inside a flame-tube burner using a special 
water-flushed probe. The soot is measured at a series of points at each burner and 
upon occasion gaseous constituents NO, CO, hydrocarbons, etc. have also measured. 
Four geometries of flame-tube burners have been studied, as well as a variety of 
different fuels. The main results of the report can be put into categories. First, we have 
studied the role of upstream geometry on the downstream pollutant formation. We have 
found that the amount of soot formed is particularly sensitive to how areodynamically 
clean the configuration of the burner is upstream of the injector swirl vanes. Second, 
we have studied the effect of pressure on soot formation. We have chosen to interpret 
the main effect of changing pressure as a variation in the Reynolds number of the flow. 
We have found that beyond a certain Reynolds number, the peak amount of soot 
formed in the burner is constant. Third, we have studied a variety of different fuels in 
the burner to determine a relationship between the properties of the fuel and the 
amount of soot formed. We have found that the hydrogen to carbon ratio of the fuel is 
a rough indicator of the peak amount of soot formed in the burner — although one can 
distinguish, in our data, the effects of different fuels having the same hydrogen to 
carbon ratio. That is, some fuels produce more soot than others evtn when the 
hydrogen to carbon ratio is the same. Finally, a discussion of these results in the line of 
current literature is provided with tentative interpretation of the data. 


I. 


INTRODUCTION 


This research has been motivated by a variety of interests. It was recognized that 
very little firm information is available on the kinetic scheme of soot formation in a 
typical gas turbine combustor, so that theoretical models of the process, which would 
necessarily involve uncertainties in turbulence estimates, kinetics, local stoichiometry, 
and local temperature of the flow, seemed unpromising. Instead, this report discusses a 
series of experiments, which have arisen as we have more carefully defined the 
conditions of a laboratory system so that it behaves in a similar fashion to the 
primary zone of a combustor in a large gas turbine. 

The geometry of these experiments has been fixed downstream of the injector. 

It consists of a fuel injector (air assist, or pressure type atomizer) on the centerline of a 
circular cross section flame tube. The atomizer is supported by an axisymmetric set of 
45° swirl vanes. 

Initial efforts were made to develop a reliable sampling system, to measure the soot 

( 2 ) 

concentration and polycyclic hydrocarbons (PCAH) in the flow. 

Even though the geometry of the burner seemed straightforward, we soon found out 
this was not the case. The axisymmetric geometry did not produce an axisymmetric 
mean concentration pattern for soot, hydrocarbons, CO, etc. A variety of 
geometric changes were studied; In each soot concentrations were measured. We have 
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shown that the geometry of the hardware upstream of the burner needs careful design 
If axisymmetric, uniform concentrations of soot are produced. Section II of this report 
details the geometries studied. 

The results of the purely geometric studies is presented in Section III. In that section we 
discuss how changes upstream of the burner make substantial changes in the emission 
generated. 

Due to the atmospheric pressure and low velocities used in our first burners, the 
Reynolds numbers were low, about 5,000, based on hot flow exit conditions. The 
turbulent properties of the flow could be expected to change with increasing density 
(pressure), velocity, and size - all of which increase tl.c Reynolds numbers. In Section 
IV scaling laws for soot formation are presented. These laws are as yet tentative. 

In addition. Section IV deals with a variety of results: The location of PCAH compounds 
in the primary zone; the sooting behavior of benzene - heptane mixtures; the effects of 
varying stoichiometry on soot formations; and the sooting effects of varying 
hydrogen/carbon ratios in blends of jet engine fuels. 

The long range goal of this research is to be able to reliably simulate, in a laboratory 
burner, the chemical and fluid processes which occur in a gas turbine combustor. To 
do this, we have developed preliminary scaling laws discussed in the last section 
(Section V) which suggest how soot levels vary with fuel and stoichiometery. In the 
long run, this engineering approach should lead to a Defter understanding of how to 
design and adapt combustors to new and changing fuels. 
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II. 


GEOMETRIES TESTED 


There are four geometries which will be discussed in this study, which we shall label A, 
B, C and D. Figure 1 shows the geometries A and B. In geometry A, the burner is 
tubular, 53 cm long, combustion air is supplied at atmospheric pressure, and is passed 
through k5° stationary swirl vanes. The burner diameter is 12.7 cm. Fuel is supplied 
through an air assist atomizer (DeLavan n.odel 3060-9-11, see Fig. 2) on the burner 
axis. The air supply for the atomizer is independent of the main air supply. Perhaps 5% 
of the total air flow in the burner passes through the atomizer. The atomizing pressures 
vary from 18^ KPa (12 psig) to 239 KPa (20 psig). Ignition is accomplished with a pair 
of standard oil burner electrodes. 

Geometry B is identical to geometry A, except the diameter is 7.6 cm, hence allowing 
for higher flow velocities, while using the same compressor. 

The upstream configuration of geometries A & B is shown in Figure 1. Immediately 
upstream of the swirl vanes, the pipe diameter is the same as the burner diameter. 
Within one diameter of the swirl vanes, the fuel and air are fed into an air blast injector 
through standard plumbing hardware. In the plane of the air and fuel fittings, the 
blockage in the burner is about 30% of the area of the tube, for the 7.6 cm diameter 
burner. 

Geometry C is the same as Geometry B downstream of, and including, the swirl vanes 
and atomizer. However, upstream of the swirl vanes, the pipe diameter Is 15.2 cm, so 
the blockage due to fittings is about ^/k that of Geometry C. In addition, the main air 
flow passes through a large settling chamber before entering the pipe section. See 
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Figure 3. This system can be run premixed by evaporating fuel in the settling chamber 
at a known rate, as well as operating in the standard heterogeneous mode. 

Geometry D is a 7.6 cm diameter tube burner, with 45° swirl vanes, using a pressure 
Injector. The fittings, upstream of the swirl vanes, used to carry fuel to the nozzle are 
specially designed and fabricated, and consist of aerodynamically streamlined struts 
positioned some 3 diameters upstream of the nozzle. Fuel flows through the struts to a 
tube, located on the centerline, and hence to the nozzle. Geometry D is capable of 
being pressurized to about 7 atmospheres. The mass flow rate is controlled by a sonic 
orifice downstream of the flame tube. See Figure 4. 

Gas and soot samples In each of these geometries are collected using a unique, water 
flushed probe. Figure 5 shows this probe, which, when it is used to collect soot, is 
flushed by part of the cooling water flowing down the center tube. The probe flush 
water is turned off when gas samples are taken, but small amounts of leakage prevent 
reliable measurements of water vapor to be made. 

The soot collected Is deposited on a fiberglass filter, which Is desiccated and then 
weighed to determine the amounts of soot collected. As the sampling flow rate is 
known, the resulting figure can be translated into soot mass concentrations Inside the 
burner. 

The gas concentrations measured, CO, CO^f NO^ HC and are measured by 
commercially available instruments mounted on a gas cart. 
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GEOMETRIC EFFECTS 


The first indication of the complexity of the How came from measurements of the 
contours of soot in geometry B. Typical examples are shown in Figure 6, which gives 
the estimated soot contours 2 diameters downstream of the swirl vanes and injector. It 
is seen that the soot pattern is very far from axisymmetric. For this soot pattern, the 
injector was an air assist atomizer using Kerosene run with 15 psig atomizing pressure, 
and fuel equivalence ratio of 1.0, Study of the gradients in soot closer to the burner 
produced the data shown in Figure 7 taken along one radius, (tl fixed in cylindrical 
polar coordinates) for Z/D = 1.33, 

Since we wished to determine if this effect was repeatable, without the long data 
collection time for soot measurements ( 3 hours/data point), we measured gas 

concentrations (CO, CO^# ^2 hydrocarbons). An informative way to plot the 
data is to use the water gas equilibrium constant (Fig. 8) to reduce the measured 
concentrations to values of r . The results of repeated runs is shown in Figure 9 at 
Z/D of 1.33. Although there is about a 15% scatter in the data, the effect is clearly 
reproducible. Apparently, the rich area of the flow is formed near the injector. But the 
effects of the nonuniformity are reflected In nonuniform soot concentrations some 
distance downstream. 

To define this effect further, the upstream end of geometry B was rebuilt, as described 
In the previous section. In this geometry C, we compared the maximum centerline 
value of soot concentration (at about Z/D = 1.5) for the B and C geometries. This is 
shown in Figure 10 for benzene. An additional feature of the geometry C is the ability 
to run premixed, this condition is also shown for heptane in Figure 1 1. 


Keeping In mind that the geometry and flows through the Injector and swirl vanes are 

Identical In the old geometry (B) and the new geometry (C) In Figure 10, we can see 

that substantial differences In soot level can be caused by different upstream 
(5 3) 

geometries. ’ The Injector remained In place during the premixed tests. The fact 
that premixing reduced the soot levels to zero over the range 0.4 <<|><1.6 Is consistent 
with thermodynamic equilibrium calculations which suggest that no solid carbon should 
exist at these values of I*, and that soot formation Is related, In some way, to the non- 
uniform mixing of fuel and air. 

Turning now to geometry D, the nonuniformity In soot was controlled by careful 
attention to aerodynamic detail upstream of the nozzle and swirl vanes. Figures 12, 13 
and 14 show measured values of soot for various positions. These data Indicate that 
the soot and pattern Is axisymmetric if both the upstream and downstream geometry 
Is axisymmetric. 

In geometry D we made a special study of the role of the blockage such as could be 
caused by an upstream obstruction. In the aerodynamically clean geometry which gave 
the results shown In Figures 12, 13 and 14, we deliberately changed the upstream 
geometry by placing a blockage plate over one of the swirl vanes (shown schematically 
In Figs. 15 and 16). 

In performing tests to determine the role of this blockage, a limited test grid was set 
up, comprising of a measurement along the center line of the burner and measurements 
at four azimuthal angles, 0= 0®, 90®, 180® and 270®, at a radius of 2.5 cm. The 
sampling grid went from iVj diameters downstream of the injector to 2, 2 % to 3, 4, 5, 6 
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and 7 diamaters downstream. A large number of data points were generated. To cut 
down on the collection time, measurements of the value of 4> were determined through 
exhaust gas sampling - which Is much faster than soot collection. 

It's convenient in summarizing this data to extract from the data the data points which 
show the effect of the blockage from those which do not. The data points which show 
the effect of the blockage are called anomolous points. The figures which we show 
next separate the anomolous data points from the other data points. That is, only a few 
measured values in each cross section look unusual when only one vane of the swirl is 
blocked. All the other points are called normal. For example, in Figure 17 we show 4> 
at a radius of 2.5 cm, for various Z/0, the normal data points as a function of axial 
location both with the blockage plate in and the blockage plate out. The value of the 
first rises and then levels off at about Z/0 = 4. The final value, = .8, for Z/D greater 
than 4, is the overall stoichiometry of the burner. In the next figure. Figure 18, the 
same data is shown for the center line. The center line starts veiy rich, because at 
Z/D = 1, the fuel spray is very rich on the burner axis. Again by Z/D - 4, the flow 
has become well mixed to the extent that the center line value of is essentially the 
same as that of the inflow fuel and air. [These values of are determined using the 
water gas reaction in the measurements of gas constituents.] 

Excluding the anomolous points, Figure 19 shows the oxygen concentration as a 
function of distance downstream and one can see that for Z/D > 4 the equilibrium 
(corresponding to<^ of 0.8) oxygen concentration of about 4% is obtained. 

We now can ask how far downstream the nonuniformify due to blockage persists. This 
has to do with the anomalous data points. For example. In the next figure. Figure 20, 
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we show the value o* 4) for various azimuthal positions of the probe location at Z/0 ■ 5. 
Although the error bars overlap at e ■ 90° • 160°, and 270* they do not at @ > 0*. The 
error bars are from a number of reported points. From this we can conclude that <tt 5 
diameters there Is an effect of blockage which can be seen in the contours of 
stoichiometry. A similar result occurs for all diameters less than 5. The final figure in 
this section. Figure 21, shows that at 6 diameters, all the data points essentially 
overlap. Other data show the same thing for Z/D of 7. Although the grid Is not fine 
enough to determine whether the nonuniformity rotates with the supposed hot flow swirl 
velocity set up by the swirl vanes, the data do support the conjecture that the 
nonuniformity persists a substantial distance, 5 diameters, downstream. \Mthin the 
accuracy of the data, it is between H and 5 diameters downstream that the flow 
becomes well mixed. 

These data support the idea that a nonuniformity upstream can generate a defect or a 
change in stoichiometry which persists until combustion is essentially complete. 
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IV. 


CHEMICAL EFFECTS 

The firs! measurements of chemical constituents in the turbulent, swirling flow 

characterized the soot formation process in the geometries which we presently label A 

( 2 ) 

and B. These results were reported by Prado, et al In 1976. The essential result will 

only be summarized here. Figures 22 and 23 show the soot concentration of 

geometries A and B respectively, (Figures 2 and 3 of Prado, et al). Compare these 

( 2 ) 

results with Figure 2k showing the concentration of PCAH (Figure k, ). Polycyclic 
hydrocarbons, (PCAH), were found closer to the nozzle or Injector than the peak in 
soot concentration. In analyzing the structure of PCAH some 30 cyclic compounds 
were identified, with the molecular weight increasing with distance from the nozzle. In 
geometries A and 6, where Prado, et al., got their results, the flow is nonuniform due to 
the geometry upstream of the swirl vanes, and the Reynolds numbers are small, an 
effect which will discussed In a moment. 

(8) 

The next study of chemical effects was done In geometry C. The peak amount of 
soot formed on the axis is shown in Figure 25 for a variety of benzene - heptane 
mixtures. The point shown on this chart, "Kerosene”, is so iocated that its C to H ratio 
matches that of the appropriate benzene - heptane mixture. One can see that there is 
substantial difference between the peak amount of soot formed from benzene - 
heptane mixture and from Kerosene, at the same C to H ratio. The Reynolds number of 
the flow is low, as this burner operated at atmospheric pressure. 

Another chemical effect which is important is that of stoichiometry. Our current 
understanding of the role of the overall stoichiometry of the burner runs something like 
this. Only a fraction of the eddies in the (low form soot. In particular, if the burner is 
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running with an overall equivalence ratio of near 1, there will be some probability 

distribution of for each eddy, around the mean value of the stoichiometry. From 

thermodynamic calculations, we believe that only eddies which are richer than a 

stoichiometry of 2 will form soot. Hence, only a small fraction of the eddies In the flow 

are rich enough. However, as the overall value of 4 > is varied toward the richer side 

there will be a larger and larger fraction of eddies rich enough to form soot. This effect 

( 6 ) 

shown In Figure 26 for geometry B. It is interesting and Instructive to compare this 

(9) 

result with that of MacFarlane which is shown on the same curve. MacFarlane's 
data was taken at 6 atmospheres whereas the data shown In Figure 26 was taken at 
atmospheric pressure. One can see that there Is not a major difference between the 
two results even though MacFarlane's geometry was not the same as geometry B. 

We now turn to the results of chemical effects which we obtained in geometry D, 
the pressurized burner. Our first efforts were directed to the following experiment: 
The burner was run at a fixed overall stoichiometry of 1.25 with Kerosene as a fuel. 
The upstream and downstream geometry of the burner was held fixed. In this system 
the pressure of the burner was changed systematically by changing the size of the 
sonic orifice downstream. In this way the Reynolds number of the flow was changed 
due to the changing density of the flow. We then ask how the peak value of soot 
varies with pressure or Reynolds number. These results are shown In Figure 27. The 
middle two points In the chart are for the same Injector but the data point at a Reynold' 
number of approximately 12,000, and the point at the Reynolds number approximately 
3,500 are for two other injectors. However, this initial result seems to suggest that 
above a certain Reynolds number, the peak amount of soot formed was Independent of 
Reynolds number. Each plotted point is the average of about 3-6 measurements. The 
accuracy is - 10%. 
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To remove the effect of changing injectors, the same experiment was repeated using 
propane flowing out through a jet coaxial with the center line of the burner in the same 
location as the liquid fuel injector. In this experiment there was no change in tne 
Reynolds numbers, swirl numbers, or geometry over the range of pressures studied. 

The results are shown in Figure 28 and it is clear that there is a leveling off at a certain 

Reynolds number on the order of 15,000. These results imply an absence of pressure 

(9) 

scar j of the peak soot loading. This is contrary to the results of MacFarlane (1974), 
a point which we will take up in the discussion. 

Figure 29 and 30 and all the remaining figures in this section deal with measurements 
made at Reynolds numbers high enough to be in the range where the peak amount of 
soot is independent of Reynolds number. In Figure 29 and 30 we show four different 
fuels. Jet A, Decaline, Tetralin, and Xylene. Notice that the peak value of soot is 
always at the same location in the burner. Within the accuracy of the data, the rise 
from Z/D = 0 to the peak value of soot and the subsequent decay afterwards has 
essentially the same slope, for the different fuels. For these reasons, both the data in 
•figures 29 and 30 and also in the preceding two charts are use only the peak value of 
soot as an indicator of the total behavior of soot in the burner. 

The last chart shows the peak loading of soot scaled with percent hydrogen by weight 
for the various blends tested. It is commonly believed that the C to H ratio is a good 
way to scale the amount of soot formed for different fuels (10) and this figure, Figure 
31, supports that idea. However, one can see from the data that there is difference 
between the amount of soot formed and the character of the blend. We believe that 
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these closely related blends might be typical of what would get from broad spec fuels, 
such as the ERBS fuel. Hence, if one scales entirely by percent hydrogen one 
could Imagine the different fuels would give somewhat different levels of soot In a given 
combustor. 
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V. 


DISCUSSION OF RESULTS 


The major part of this discussion will be centered around the interpretation of the 
Reynolds number scaling presented In the previous section. However, It should be 
emphasized at the outset that these data are not complete because the variation with 
Reynolds number was achieved only with a variation in pressure. Ideally, these 
experiments need to be repeated In which a given Reynolds number was achieved by 
varying either flow rate or pressure. In the experiments which we report the pressure 
and flow rate were varied simultaneously rather than independently. 

We start by attempting to make a nondimensional argument for how the flow geometry 
and combustion process in the burner might scale. A number of the important 
nondimensional parameters are known. First of all, there is the stoichiometry 


V = 


mass flow fuel/mass flow air 


(mass flow fuel/mass flow air)/ 


stoichiometric 


( 1 ) 


Typically, in our experiments, 4' varies from 1.5 to 0.6. However, in the section on 
chemical effects, most of the s' values were near 1. 


A second scaling parameter is a Reynolds number: the ratio of viscous to inertial 

forces. It is Known from standard non^reacting flows that the Reynolds number is 
important in determining the tuibulent motions of the flow. When the Reynolds number 
of a flow Is high enough the large scale turbulent mixing processes no longer depend on 
viscousity. That is, the turbulent diffusivity for the large scale eddies is Reynolds 
number independent. For a pipe flow with no reactions, the friction on the wall does 
not vary Reynolds number for P*' 10-20.000. 

Re (2) 




In present experiments the Reynolds number ranged from a few thousand for geometries 
A, B and C up to approximately 20,000 In various circumstances for geometry D. 

Another important nondimenslonal number characterizing the flow is the swirl number. 

( 12 ) 

Tnis number has been described by Beer and is a ratio of the angular momentum 
to the axial momentum gning through the swirl vanes. Its definition Is (see Beer, Sec. 
1 . 1 ) 



Un2 rdr 



l-puairdr 


(3) 


A typical value for the swirl number in all the experiments reported in this report is 0.8. 
Varying the swirl number changes the amount of recirculation just downstream of the 
injector swirler geometry. 


Fuel properties are more difficult to scale. First of all, there may be many reactions for 
different parts of the combustion process and, secondly, different fuels may have 
different reaction schemes. However, the hydrogen to carbon ratio has been suggested 
by the results of Section IV as a parameter. It is simply the percentage weight of 
hydrogen to the percentage weight of carbon in the fuel. We know that this scaling is 
not completely correct because different fuels which have the same C to H ratio give 
slighiiy different results. However, the early results shown in Figure 25 on benzene - 
heptane mixtures In geometry C should not be interpreted as a k ck of hydrogen to 
carbon scaling of soot loading. Those first experiments were done at too low a Reynolds 
number. For the want of a better scheme, we shall assume that X scales the fuel 
properties. 
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V _ mass Hydrogen , % 

“ mass Carbon 

Finally, we can make an estimate of the ratio of the flow time divided by the reaction 
time to form soot. To estimate the reaction time, we look up low pressure data for 
burners which has been summarized by Prado These data show a time of about 1 
millisecond to reach a peak soot concentration in a laminar flow burner at 20 tor. We 
shall assume that the soot reaction time scales inversely with pressure, as does the 
flame thickness. This is based on the fact that the soot seems to be formed at the 
same relative position in the laminar flame, independent of pressure (1). Using this 
estimate, we can calculate the nondimensionat ratio of the flow time to the reaction 
time. 


flow time 
reaction time 


~ •» 760 mmHg * 5 atm P 3 . . 

10~ sec. (20)mmHg U 


This calculation indicates the flow time or turbulent mixing time is much longer than the 
reaction time. This supports the idea that the turbulent mixing process controls the 
amount of soot formed. Hence that there should be a Reynolds number scaling of the 
soot, as the mixing process scales with Re. 


However, we have no definitive argument about why the peak values of soot should or 

(9) 

should not scale with pressure. MacFarlane's data, which suggests a strong 
pressure dependence of soot amount, is somewhat deficient in this regard because, in 
the same paper, he showed that the flow geometry of his burner system changes with 
pressure. As we have found, (Sec. Ill) small changes iri geometry can make major 
changes in the amount of soot formed. However, MacFarlane's data can be 
contrasted with the data of Kadota et al for droplet burning, at various pressures. 
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They show that when the ambient pressure of the atmosphere around the drop rises to 
the point where the burning wake of the droplet becomes turbulent, the mass of soot 
produced from burning a droplet essentially remains constant as the pressure Is further 
increased. Although the present work does not refer to droplet burning, a similar 
scaling effect may occur In our fully turbulent flow. 


If it turns out that future studies confirm the scaling suggested here, interesting 
possibilities for the practical application of this similitude, based on the five 
nondimensional numbers listed, seem important. For example, it might be possible, by 
setting up an appropriately scaled swirler and injector, to simulate the amount of soot 
formed one or two diameters downstream of the swirler in a full scale engine. 
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Calculation of equivalence ratio, 0, from 
Values of [NO>], [CO], [COa], (HC], and [Oa] 

Equation of combustion 

0C>HyOi + 4.76(Air)-»>AiCO + AaCOa AaHiO 

+ AaHi + AlOa + A»CfH» 
+ Ar(NO«) 

Measured: Dry concentrations of CO, COa, Oa, & NOa 
Wet concentration of HC 

Goat: Calculate equivalence Ratio, 0 

Method and Assumptions: 

—Water-gas equilibrium constant 
j^_AiAi_3 8 
AaAa 

—Calculate [HaO] using [CO] & [COa] and K. 

—Wet concentrations found from dry concentrations 
[ ]w = [ ]o[1-[HaO]] 

— ForHC x_r 

7 7 


Fig. 8 Method of Determining Equivalence Ratio of Sampled Gas as 

a Function of Species Concentrations. 
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and Without Blockage, Configuration 0. 
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